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ABSTRACT 

Many competing linear instabilities are likely to occur in astrophysical settings, and it is important 
to assess which one grows faster for a given situation. An analytical model including the main 
beam plasma instabilities is developed. The full 3D dielectric tensor is thus explained for a cold 
relativistic electron beam passing through a cold plasma, accounting for a guiding magnetic field, a 
return electronic current and moving protons. Considering any orientations of the wave vector allows to 
retrieve the most unstable mode for any parameters set. An unified description of the Filamentation 
(Weibel), Two-Stream, Buneman, Bell instabilities (and more) is thus provided, allowing for the 
. exact determination of their hierarchy in terms of the system parameters. For relevance to both 

■ real situations and PIC simulations, the electron-to-proton mass ratio is treated as a parameter, and 
numerical calculations are conducted with two different values, namely 1/1836 and 1/100. In the 
system parameters phase space, the shape of the domains governed by each kind of instability is far 
from being trivial. For low density beams, the ultra-magnetized regime tends to be governed by either 
the Two-Stream or the Buneman instabilities. For beam densities equalling the plasma one, up to 
four kinds of modes are likely to play a role, depending of the beam Lorentz factor. In some regions 
of the system parameters phase space, the dominant mode may vary with the electron-to-proton mass 
ratio. Application is made to Solar Flares, Intergalactic Streams and Relativistic shocks physics. 

O ' Subject headings: instabilities cosmic rays gamma rays: bursts - shock waves 

. 1. INTRODUCTION 

Weibel, Filamentation, Two-Stream, Bell or Buneman instabilities are ubiquitous in astrophysics. They are involved 
CN| ' in the physics of Solar Flares where relativistic electron beams are assumed to lose their energy through beam-plasma 
instabilities, producing hard X-ray emissions (Karlicky et al. 2008; Karlicky 2009). Filamentation, or Weibel, insta- 
bilities could also be responsible for the birth of cosmological magnetic fields (Schlickeiser & Shukla 2003; Schlickeiser 
1 2005; Lazar 2008; Lazar et al. 2009) as unstable particle streams through the intergalactic medium can magnetize 
an initially un-magnetized system. Such instabilities could also play an important role in explaining the origin of a 
variety of high energy photons sources including Supernova Remnants, Active Galactic Nuclei, Gamma Ray Bursts 
Cn or Pulsar Wind Nebulae (Piran 1999; Gedalin et al. 2002; Piran 2004; Waxman 2006; Stockem et al. 2007). Within 
| some of these systems, it is assumed that cosmic rays are accelerated through shocks (relativistic or not) while the 
ON . instability generated upstream by their interaction with the interstellar medium provides the magnetic turbulence 

■ eventually responsible for synchrotron radiation emissions (Medvedev & Loeb 1999; Silva et al. 2003; Nishikawa et al. 
^ I 2006; Frederiksen et al. 2004; Milosavljevic & Nakar 2006; Lemoine et al. 2006; Niemiec et al. 2008). 

Regardless of the context, the typical structure investigated consists in a beam-plasma system initially both charge 
and current neutralized. Initial charge neutralization implies the inclusion of positive and negative species, while 
[ current neutralization demands at least two streaming species. Streams are often considered as opposed, but they can 
be parallel when streaming species are of opposite signs. For example, a pair beam does not need anymore beams to be 
current neutral. Finally, accounting for an external magnetic field Bo allows the system to be relevant to a wide class 
of astrophysical problem. The simplest case consists in a Bo parallel to the flow(s) but normal or oblique orientations 
have also been considered (Fraix-Burnet & Pelletier 1991; Dieckmann et al. 2008; Bret & Dieckmann 2008a; Sironi & 
Spitkovsky 2009). 

We consider here the generic system formed by a cold electron beam of density rib streaming at initial velocity v;,, 
with Lorentz factor jb = (1 — v 2 /c 2 ) -1 / 2 , over a cold electron/proton plasma. The beam current is neutralized by an 
electronic return current of density n p and velocity v p such that n^Vf, = n p -v p . Ions are initially at rest with density 
rii = rib + %>■ Finally, the model includes a flow-aligned magnetic field Bo || v&. Even for such a simple system, linear 
stability analysis is intricate because unstable modes are numerous. The intent of this paper is to clarify this issue and 
determine the fastest growing mode for any given set of parameters, by implementing an exact model encompassing 
every basic instability. In order to keep the results tractable, emphasis is made on the simplest possible model with the 
lowest number of free parameters. As explained in the sequel, these requirements demand a beam-plasma system with 
mobile ions and a guiding magnetic field, while kinetic effects will not be investigated here. A cold fluid model is thus 
exactly solved, rendering every possible coupling between unstable modes. The dispersion equation arising from the 
exact dielectric tensor is analyzed without calling on the electrostatic (k || E), or the purely electromagnetic (k _L E), 
approximations. Such kind of calculation is mandatory if the Two-Stream or Buneman electrostatic modes are to be 
described within the very same formalism than the electromagnetic Filamentation or Bell-like instabilities. 
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Let us list the possible unstable modes by progressively "assembling" the system. The first block here is the 
un-magnctized background plasma with fixed protons. Because we take it cold, it is stable against every kind of 
perturbations. Note that some temperature anisotropy would make it Weibel unstable (Weibel 1959), the fastest 
growing modes being found for wave vectors perpendicular to the high temperature axis (Kalman et al. 1968). Kinetic 
effects can thus drive this first component unstable even before we add anything else. 

Let us now "add" the electron beam and let the plasma electrons establish the return-current (protons are still fixed) . 
To simplify the discussion, we consider here a diluted beam with rib <C n p (the rest of the paper also deals with higher 
beam densities). As is known, the resulting beam-plasma system is unstable. Perturbations with wave vector k || v& 
arc unstable for < k < uo p /v b and define the so-called Two-Stream Instability (Bohm & Gross 1949). Perturbations 
having k _L v& are also unstable for any k and pertain to the Filamcntation instability 1 . Finally, perturbations with 
k neither parallel nor perpendicular to the flow are also unstable (Famberg et al. 1970; Califano et al. 1997, 1998a, b) 
and may even be the fastest growing modes in the diluted relativistic beam regime (Bret & Deutsch 2005; Bret et al. 
2008). 

Let us now give the protons the possibility to move. Every aforementioned unstable modes are still unstable as the 
two counter-streaming electron beams still interact. But both electron currents can now interact with the protons, 
giving rise to Buncman unstable modes (Buneman 1959). A closer look at the situation shows that Buneman modes 
arising from the beam/protons interaction merge with the beam/return-current modes (Bret & Dicckmann 2008b). 
But the reiwri-current/proton interaction results in unstable Buneman modes reaching their maximum growth rate 
for k ~ (n p /nb)Ldp/vi, and k || v&. Within the diluted beam regime, rib -C n p implies unstable modes with much shorter 
wave length than the Two-Stream modes. Indeed, it has been found than the Buneman modes can compete with the 
Two- Stream ones if the beam is relativistic enough. 

We finally make our system complete by adding a flow-aligned magnetic field. Here again, the new ingredient 
brings in more unstable modes without necessarily stabilizing the previous ones. To start with, electronic Cyclotron 
and Upper-Hybrid modes are destabilized, adding new branches to the unstable spectrum for any orientations of k 
(Godfrey et al. 1975). Also, Alfven modes resulting from the combination of the vibrating protons in the magnetic 
field can also be destabilized. This unstable modes where first pointed out by Bell by means of a MHD formalism 
(Bell 2004, 2005), and their description was later extended through the kinetic one (Zweibel 2002; Reville et al. 2006). 

Our simple magnetized beam/plasma system with moving protons in eventually Two-Stream, Filamcntation, 
Oblique, Cyclotron, Upper-Hybrid, Buneman, Bell. . . unstable! As the system is released from equilibrium, every 
instabilities are triggered, and the outcome of the linear phase is mostly shaped by the growth of the fastest growing 
one. This most unstable mode can be found for an oblique wave vector so that the search of the dominant mode re- 
quires the implementation of a model capable of describing any unstable mode for any orientation of the wave vector. 
Note that recent works involving ion beams also evidenced such kind of modes (Niemiec et al. 2008; Ohira & Takahara 
2008). We thus now proceed to the elaboration of the simplest possible model incorporating all the aforementioned 
instabilities. 

The unstable modes described here are just too numerous to be detailed one by one. On the other hand, each one 
needs to be discussed since the present aim is precisely to determine their hierarchy The full wave vector dependance 
of each mode is therefore skipped, the focus being set on the most unstable wave vector with the corresponding growth 
rate. Additionally, modes which have not been found to govern the system for any given sets of parameters are only 
briefly discussed in order to keep the presentation tractable. Still for clarity, the names of the numerous modes have 
not been abbreviated. Finally, analytic expressions in terms of the magnetic field parameter fig defined by Eqs. (4) 
are essentially derived for Qb > 1, while results are also presented for the un-magnetized regime. 

The paper is structured as follow: the analytic model is explained in the next section. Unstable modes found for 
flow-aligned wave vectors are listed in Section 3 before an overview of the 2D unstable spectrum is given in Section 4. 
The key results of the article are exposed in Section 5 where the fastest growing mode is determined in terms of the 
system parameters. In this respect, Figures 5 and 8 can be considered as the main results of the papers. They show 
which kind of mode governs the linear phase of the beam-plasma system in terms of the parameters. Application is 
then made to various astrophysical settings in Section 6, before the final discussion and conclusion. 

2. ANALYTIC MODEL 

The model relies on the cold relativistic fluid equations for the three species involved. The calculation follows the 
lines of previous ones (Califano et al. 1998b; Kazimura et al. 1998) except that protons are allowed to move while a 
static flow-aligned magnetic field is accounted for. The basic equations include Maxwell's equations and, 

^+V.(n 3 v J ),0, (1) 



dt 



E | Vj - x (B + B ) 



(2) 



where j — b for the electron beam, j — p for the plasma return-current and j = i for plasma ions (protons). Here, 
Pj = "fjtnjVj and qj are the momentum and electric charge of specie j. 

1 The Filamentation instability is singled out from the Weibel one. In its original context, the Weibel instability results from a temperature 
anisotropy with no drift. Filamcntation instability involves relative streaming of various species. In the present context, both instabilities 
maybe be disconnected as one Filamcntation stable, or unstable, beam can interact with a Weibel stable, or unstable, plasma (Bret et al. 
2004, 2005). They can also interfere with each other (Bret & Deutsch 2006; Lazar et al. 2006; Lazar 2008; Lazar et al. 2009). 



Weibel, Two-Stream, Filamentation, Oblique, Bell, Buneman. 
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TABLE 1 

Main unstable modes in the diluted beam regime. 



Modes 


Z x 


z* 


Growth rate 


T X X ^T X y 







<x ^/aRftg'"" 











Two- Stream 





1 


#« 1/3 /7 6 


Buneman 





1/q 


V% pl/3 

,4/3-"" 


Obliqucgoi f^s = 
UHL***, n B > 1 
Oblique, Qb > 1 
UHL***, Q B > 1 


CO 
CO 

^/i + n| 

CO 


1 

^/i + n| 


AWti) 1/3 

K /3 /76 



For a given , the growth rate is a function of Z z . The table only mentions the Z z and the growth rate pertaining to the fastest growing 
mode. Magnetized mode are reported for Qg > 1. The transition between Qg = and 1 is described in Godfrey et al. (1975) for R = 0. 
"Bounded by Eq. (13). 
***UHL stands for "Upper-Hybrid-Like". 

Although lengthy, the derivation of the dielectric tensor is quite standard. The conservation, Maxwell's and Euler's 
equations are first linearized assuming every quantity slightly departs from equilibrium like exp(zk • r — iujt) where 
i 2 = —1. Since Bo || v p while Qj n j = 9i n i v j = 0; the charge and current neutral equilibrium state considered 
exactly fulfills the full set of equations. We can write k — (k x , 0, k z ) without loss of generality by virtue of the axial 
symmetry with respect to the beam and magnetic field axis (Godfrey et al. 1975). We choose the z axis for the 
direction of the beam and of the magnetic field, having therefore = (0,0, Vb), v p = (0,0, v p ) and B = (0,0, -Bo)- 
The linearized conservation and Euler's equations are first used to express the perturbed total current Ji in terms of 
the electromagnetic field. The first order magnetic field is then eliminated through Bi = (c/w)k x Ei and the resulting 
expression of the current is inserted into the usual combination of Maxwell-Faraday and Maxwell- Ampere equations, 
namely 

c 2 4?7r 

Ux(kxE 1 )+E 1 + Ji(E!)=0. (3) 

U z U 

The electrostatic (kx Ei = 0) or purely electromagnetic approximations (k-Ei = 0) would here result in a simplification 
of the k x (k x Ei) term. Although such options may be useful when focusing on one give mode, the risk at this stage 
would be to loose informations about possible dominant modes which would not fit in. A Mathematica Notebook 
has been implemented to symbolically compute the tensor (Bret 2007). It is expressed in terms of the dimensionless 
variables, 

w ku b n b v b u b m e 

x= — , Z = , a=— , 13=—, Q B = — , R = — , (4) 

uj p u p n p c Up n%i 

where u 2 — 4Trn p q 2 /m e is the electronic background plasma frequency, u b = \q\B /mc the non-relativistic electronic 
cyclotron frequency and m e , m, the electron and ion (proton) mass respectively. The Alfven velocity va can be 
expressed in terms of the variables above as, 

— = ^R(l + a)n B , with v A = /A B ° ■ (5) 

The dielectric tensor has the form, 

T xy T yy T yz , (6) 

1 xz 1 yz 1 ZZ J 

where z* is the complex conjugate of z. Tensor elements are reported in Appendix A in terms of the dimensionless 
variables (4). In the limit of zero magnetic field, only the square of the charges appears in the tensor expression. 
As a consequence, the results in such case also apply to two counter-streaming pair plasmas (Jaroschek et al. 2005; 
Ramirez-Ruiz et al. 2007). 

The most general expression of the dispersion equation then reads detT = and necessarily encompasses every 
possible unstable modes. As previously mentioned, we here emphasize the search of the most unstable mode \<lm of 
the unstable k spectrum. The function kjf can only depend on the bcam-to-plasma density ratio a, the beam Lorentz 
factor 7b, the magnetic field strength parameter VLb and the mass ratio R. Once the former has been fixed to 1/1836, 
or any higher value more suited to comparison with PIC simulations (Jaroschek et al. 2004; Spitkovsky 2008a; Karlicky 
2009), the most unstable mode can only be a function of (a, 7&, f2#) e [0,1] x [l,oo] x [0, oo]. Note that this study is 
not limited to the diluted beam case so that the return-current flowing at v p — av b can also reach relativistic velocities. 
Tables 1 and 2 summarize the results of the two next sections as they display the main unstable modes in the diluted 
(a « 1) and symmetric (a = 1) regimes respectively. Modes properties in the magnetized regime are given only for 
Qb > 1- A detailed study (at R = 0) of the transition from Ub = can be found in Godfrey et al. (1975). 
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TABLE 2 

Main unstable modes in the high density beam regime a ■ 



Modes 


Z x 


z z 


Growth rate 











TSB*, jtR 1 / 3 < 1/2 







l/27 t 3/2 


TSB*, 76^/3 > 1/2 





Vm 


27/6 n lib 


Oblique, Qg > 1 




Q B /2 lb 




UHL** 1& 2, n B > 1 


oo 


~ ^s/7fc 


i/n B 



TSB stands for the merged Two- Stream/ Buncman mode 
UHL stands for "Upper-Hybrid-Like". 



3. FLOW-ALIGNED WAVE VECTORS 



Wc here detail the case of flow-aligned wave vector instabilities before analyzing the most general situation. Even if 
these modes do not always govern the system, such analysis is useful to get familiar with the competing instabilities, 
and make the bridge with previous investigations. When considering k = (0, 0, k z ), the tensor (6) simplifies to, 



T 

± XX 


T 

J-xy 





—T 


T 
± yy 











T 

J- zz 



(7) 



so that the dispersion equation simply reads, 

Tz Z {T xx -\- iT X y)(T xx °iT X y) 0. (8) 

This equation already contains most of the unstable modes already discussed. The T zz generates electrostatic Two- 
Stream and Buncman modes while the others are responsible for the so-called "Bell" unstable electromagnetic modes. 



3.1. Electrostatic modes 



The T zz terms reads, 



R(l + a) 



Or- 



Z z )^ b (x + aZ z f^- 

The proton contribution stems from the term cx R/x 2 , the beam one from the second term with numerator a, and 
the return-current with j p = (1 — v 2 /c 2 )~ x / 2 accounts for the last factor. The dispersion equation for these modes 
is independent of the magnetic field since the resulting electrostatic modes have k || E and only generate velocity 
perturbations along the magnetic field. This equation has been analyzed elsewhere (Bret & Dieckmann 2008b). For a 
diluted beam, the interaction of the two electron beams produces Two-Stream unstable modes with maximum growth 
rate (Mikhailovskii 1974) 



o = — 7777 lor 

2 4 /3 lh 



Z 7 = 1, 



(10) 



and the interaction of the return-current with the proton yields unstable Buneman modes at much smaller wavelength 
with, 



^|*V3 for Z z = 1/a. 



(11) 



These expressions change when the density ratio approaches unity because the beam and the return-current become 
symmetric for a = 1 (see Sec. 3.3). 

3.2. Electromagnetic modes 
The unstable electromagnetic modes arise here from the T xx ± iT xy terms of Eq. (8) which read, 

Z 2 a(Z z -x) x + aZ z Rx(l + a) 



TxX ± ^ X P 2 * l b {x - Z z ) T n B y p (x + aZ z ) T n B 



x ±Rn B ' ^ 

A straightforward calculation from the general dispersion relation for circularly polarized wave with such wave vectors 
(Ichimaru 1973; Achterbcrg 1983; Amato & Blasi 2009) yields the very same result. These modes are the present 
configuration analogs to the ones explained by Bell (Bell 2004, 2005). The main difference with the scenarios in which 
they are usually involved is that they are presently driven by an electronic current rather than by a baryonic one. 

Both dispersions equations yield here unstable modes. The most relevant feature with respect to the mode hierarchy 
is that growth rates are bounded by Ril B , regardless of the parameters involved. This does not imply the growth rate 
can go to infinity because unstable solutions exist only in a limited region of the phase space (a, jb, ^s, R)- 

The T xx — %T xy term can be partially assessed analytically neglecting the x 2 and the return-current factor. The most 
unstable modes are found at Z z ~ il B /7b, with x <~ Ril B . The maximum growth rate varies like y/a, and is always 
bounded by 

S < Rn B -n 2 B R 3 / 2 . (13) 



Weibel, Two-Stream, Filamentation, Oblique, Bell, Buneman. . . 5 



5 




0.01 0.1 1 10 



Fig. 1. — Growth rate S (u) p units) of the four unstable modes for flow-aligned wave vectors in terms of Z z = k«Vb/ui p . Parameters are 
a = 0.1, 7;, = 20, (!g = 1 and R = 1/1836. The two electromagnetic modes yield the thin plain and dashed curves from Eq. (12) with "-" 
and "+" respectively. The bold dashed pertains to the Two-Stream modes and the Buneman modes gives the plain bold one. 




Fig. 2. — Growth rate 8 (ui p units) for flow-aligned wave vectors in terms of Z z = kuvi/uip. Parameters are a = 1 for Qb = 10 and 60, 
with R = 1/1836 and 7(, = 10 3 . Here, electromagnetic modes (thin lines) can overcome the electrostatic ones (bold line). 

The T xx + iT xy factor is more involved and has been partially treated numerically. Unstable modes are here found at 
Z z ~ \a/Q.B, with maximum growth rate varying like a, and still bounded by Eq. (13). Note that these modes are 
stable for 75 < a£l 2 B so that instability in the diluted beam regime demands a low Lorentz factor. 

Figure 1 displays the growth rate of these four unstable modes in terms of Z z for some parameters triggering them 
all. One can observe how the relevant spectrum extends over 3 orders of magnitude. This point, and its consequences, 
is discussed in the conclusion. 



3.3. Non-diluted beam regime 

Because we study small perturbations of the system composed by both the beam and the plasma, the present theory 
is perfectly valid up to a = 1. However, the physical interpretation of the calculations can no longer be derived 
considering the beam is a perturbation to the isolated background plasma. The system beam+plasma has its own 
proper modes, which are quite different from the ones of the isolated plasma. Taking R — for example, the dispersion 
equation has two branches exactly given by, 

x 2 = Z 2 + ¥— g 2A, (14) 

which definitely differs from the cold isolated plasma dispersion relation x = 1, i.e u = uj p . For such symmetric 
system, the Buneman modes at Z z ~ 1/a merge with the Two-Stream modes at Z z ~ 1, and depending on 7f,i? 1 / 3 , 
the resulting mode may be governed by the electrons or the ions dynamic (see Table 2). 

For parallel wave vectors, electrostatic modes have been studied previously (Bret & Dieckmann 2008b), and two 
different regimes need to be considered depending on the product ^R 1 ^. For 7hi? 1 ^ 3 < 1/2, the fastest growing mode 
has 

S =A?2 fOT Z * = -472- ( 15 ) 

2 7b 3/2 2 7fc 3/2 
Note that this result is exact for R = 0. In the opposite limit 7 fci? 1 / 3 > 1/2, one has 

f V3i? 1/6 



2 7/ 67fc fOT Z * = ^ 2R - ( 16 ) 
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Fig. 3. — (Color online) 2D unstable spectrum in terms of the reduced wave vector Z for a = 0.1, 7(, = 2 and Qb = 3. (A) R = 0. (B) 
R = 1/100. (C) R = 1/1836. 

Regarding the electromagnetic modes, the two terms T xx ± iT xy in Eq. (8) yield the same growth rate because for 
a = 1, these two functions are equal through the transformation x — > —a;. Here again the growth rate of these modes 
is always bounded by Rfls, but instability occurs only for 75 < 1/R and < ^/Tfe. As a consequence, the most 

unstable mode for a given couple (7;,, fig) cannot grow faster than RQb < R\f% < VR- 

Figure 2 displays the growth rate of the two modes just discussed for parameter sets yielding a Bell governed system. 
Note that such situation demands a highly magnetized system. Although this plot gives the sensation that Bell-like 
modes can dominate some region of the parameters phase space, they were never found to govern the full 2D spectrum. 
This example illustrates the importance of accounting for the full unstable k spectrum instead of focusing on a given 
wave vector orientation. While Bell- like modes grow slower for non flow- aligned wave vectors, other modes grow faster 
in oblique directions. As a result, the unstable spectrum ends up governed by the latters, even if the formers can 
"reign" over the beam axis. This situation basically stems from the necessary smallness of the mass ratio R. Dealing 
with a proton beam would allow to reach R = 1, where Bell modes should definitely be found governing some portion 
of the phase space. A similar theory accounting for a proton beam would be needed to describe such situations. 

4. AN OVERVIEW OF THE 2D SPECTRUM 

Before we establish the full hierarchy of the unstable modes in terms of the parameters (a, 7^,0^), let us give an 
overview of the unstable spectrum in terms of k = (k x ,0,k z ). The basic principle here is that the unstable modes 
previously discussed for k x = remain unstable for some k x 7^ 0. The analytical analysis becomes more involved 
because all modes are coupled, as the tensor elements T xz and T yz in Eq. (6) no longer vanish. However, the list of 
unstable modes is easier to grasp in the diluted beam regime where the beam is just a perturbation to the isolated 
background plasma. We therefore single out the diluted and non-diluted beam regime in our analysis. Furthermore, 
unstable modes found in oblique direction have finite or infinite k x . Both kinds of modes can be tackled by different 
approaches, and are thus studied separately. 

4.1. Diluted beam - Z x = 00 

For the un-magnetized system with infinite ion mass, the resulting spectrum is now well understood, including within 
the framework of a full kinetic relativistic theory (Bret et al. 2008). There is a continuum of unstable modes bridging 
the Two-Stream instability (kj_ = 0) with the Filamentation instability (kn = 0). The magnetized version of the same 
system brings two additional oblique branches. For Qb > 1, the first one is found at Z z ~ f2s/7&, with a maximum 
growth rate varying like o/Qb- Still for il^ > 1, the second branch is located at Z z ~ fls/jb + i/l + with 
a maximum growth rate varying like \\faj^lB (Godfrey et al. 1975). Within the same region of the spectrum, the 
oblique modes already present at Qb = near Z z = 1 are now found at Z z ~ yT-j- & B with a maximum increment 
varying like a 1 ^ 3 /^- These later modes are commented in Sec. 4.3. At this stage, the unstable spectrum may look 
like the one pictured on Fig. 3(A) for the parameters specified in caption. 

At this junction, some comment is needed to clarify the meaning of "oblique" mode (with small "o"). Figures 3 
& 4 make it clear that unstable modes with both k x 7^ and k z 7^ are numerous. Such kind of modes have been 
so far refereed to in the literature as "Oblique" (Watson et al. 1960; Niemiec et al. 2008; Ohira & Takahara 2008), 
"electromagnetic beam-plasma instability" (Califano et al. 1998b), "coupled Two-Stream Weibel" (Jaroschek et al. 
2005) or "Mixed mode" (Frederiksen & Dieckmann 2008). The problem with such labeling is that "oblique modes" are 
here just too numerous for only one tag. In Tables 1 & 2, as in the rest of the paper, "oblique modes" bridging between 
Two-Stream and Filamentation instabilities for fig = are labeled "Obliqueso"- For Qb > 1, they evolve into what 
we presently call "Oblique" modes (capital "O"), with both k x 7^ and k z 7^ 0, and a growth rate scaling like 7 f J~ 1 . 
Finally, modes found reaching their maximum growth rate at finite Z z and Z x — 00 are labeled "Upper-Hybrid-Like" 
modes after Godfrey et al. (1975). 

The spectrum is even richer when ions are "allowed" to move, as Buneman and Bell-like modes are triggered. The 
resulting spectrum is pictured on Figs. 3(B,C) for R = 1/100 and 1/1836 respectively. An electron-to-proton mass 
ratio R = 1/100 may be relevant for PIC simulation where the mass ratio is usually incremented from its 1/1836 
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Fig. 4. — (Color online) 2D unstable spectrum in terms of the reduced wave vector Z for a = 1, 7 6 = 10 3 and H B = 60. (A) R = 0. (B) 
R = 1/100. (C) R = 1/1836. With the same parameters than Fig. 2, we see here how Bell-like modes around Z z = 0.08 do govern the Z z 
axis but not the all spectrum. Map (D), with parameters = 10, 75 = 5 and R = 1/1836, clearly evidences the oblique mode explained 
in Sec. 4.3 for Z x , Z x ~ 1. 

realistic value in order to speed up the ions dynamic, saving thus computer time. Regarding this last point, it is 
interesting to note that such a trick is possible because the growth rates associated with the moving ions increase 
with R, while the related scale length does not (see Eqs. 11). This is reflected on Figs. 3(B,C) where the dominant 
modes do not migrate with varying R. If such was not the case, some higher mass ratio could demand higher spatial 
resolution so that the computing time lost by a finer spatial resolution would not necessarily be compensated by the 
accelerated dynamic. Also on these Figures, the occurrence of the Buneman instability around Z z ~ 1/a = 10 is 
clear. The Two-Stream instability governs here the system even for R = 1/100, but the electron and ion spectrums 
are almost disconnected form each other only for R = 1/1836. Because of the smallness of the real mass ratio, both 
the time and length scales of each spectrum are disconnected from each other. In such case, the overall spectrum is 
nearly the superposition of each sub-spectrum. 

4.2. Non-diluted beam regime, a — 1 and Z x — 00 

Without moving ions (i? = 0) , the most relevant features of the unstable spectrum in this case are oblique resonances 
at Z z ~ f2s/7b and Z z ~ fls/^lb- Only the first one extends up to Z x = 00 and yields a maximum growth rate 
o~m ~ 1/^b (Bret et al. 2006). The next paragraph is devoted to an overview of the second kind of modes. This large 
Z x regime can be analytically investigated in the following way. We start deriving the dispersion equation in this limit 
setting a = 1, 7b = 7 P and Z x — > 00 in Eqs. (6). We then use the fact that unstable modes are found with 5R(w) ~ 
and develop the resulting dispersion equation to the second order in 1. In this large Z x limit, the system is found 
unstable only between, 

Z z \ = —\l 



w 2 lb 
n 2 ' 



(17) 



and 



lb ' 



(18) 



These results are exact and imply that Filamentation instability is canceled as soon as Z z \ > 0, i.e. fls > /3y/2 n /b- 

For £Ib ^> fil^ 2 , the unstable range decreases like Z z2 — Z z \ ~ (3 2 /CIb- The range of unstable wave vectors becomes 
therefore increasingly narrow, yielding a quasi-monochromatic unstable spectrum in this limit. 

When ion motion is accounted for, part of the stability domain can still be calculated exactly following the same 
guidelines. With R 7^ 0, one finds 

Z,i,r = Z zl , (19) 



and 



J z2,R 




RVL B /2 



(20) 



The maximum growth rate within this range remains close to 1/Qb, and another range of unstable wave vectors and 
similar growth rate appears at slightly larger Z z s. These new modes, arising purely from the ion motion, will not be 
investigated here although they are accounted for in the forthcoming numerical evaluation of modes hierarchy. 

Figures 4(A-C) shows the 2D spectrum for the parameters chosen for Fig. 2 with SIb — 60. The purpose of this 
choice is to observe how some oblique modes govern the full spectrum while Bell-like modes dominate the Z z axis. On 
Fig. 2, the spectrum for VLb — 60 is clearly governed by these modes at Z z ~ 0.08. The same modes, together with 
the way they evolve when leaving the Z z axis, are perfectly visible on Fig. 4(C) and one can check how they quickly 
stabilize for oblique wave vectors. 

4.3. Oblique modes at finite Z x 

Sections 4.1 & 4.2 discussed unstable modes located at Z x — 00. The cold fluid approximation typically send some 
fast growing modes to Z x = 00, whereas temperature effects tend to stabilize these small wavelengths instabilities by 



8 Bret 



QB = fiB=l nB = 3.2 £1b = 5.8 Qb = 19 




Fig. 5. — Frontiers between the domains governed by different modes in terms of a and 75, for various and R = 1/1836, 1/100. The 
arrows on the frontiers show how they evolve when fig increases. The red crosses show the parameter sets for the Solar Flares application 
explained in Sec. 6.1. 

preventing the pinching of too small filaments (Silva et al. 2002). But the occurrence of a magnetic field also triggers 
some truly locally most unstable oblique modes. The term "locally" means here that the very same mode, as defined 
by one root of the dispersion equation, reaches its maximum for a finite Z x . For more clarity, Figure 4(D) displays a 
2D spectrum with a — 1, while the other parameters have been chosen to single out this mode. 

For the diluted beam case, Table 1 mentions a maximum growth rate varying like a 1 / 3 /7b. For the symmetric case 
a = 1, Table 2 indicates a maximum oblique growth rate ~ 0.5/7^' 15 (numerical fit). Such scaling of this finite Z x mode 
bears importance consequences regarding modes hierarchy: a look at Tables 1 & 2 shows that the ultra-magnetized 
regime tends to stabilize every modes, except this one, the Buneman and the Two-Stream. While more modes compete 
for moderate the large £Ib limit is decided among these three. More details on the ultra- magnetized regime are 
given in Section 5.2. 

5. MODES HIERARCHY 

We finally turn to the determination of the modes hierarchy in terms of the parameters a, jb and Only the 
case 17 # = R = has been treated so far (Bret & Deutsch 2005; Bret et al. 2008) evidencing two main features. On 
the one hand, Oblique^ modes can govern the diluted beam regime from nb/n p ~ 0.53. For higher density ratios, 
Filamentation may dominate depending on the beam Lorentz factor 7^. On the other hand, Oblique^ modes govern 
the ultra-relativistic regime unless the system in strictly symmetric with — n p . 

After having calculated symbolically the dielectric tensor using a previously designed Mathematica Notebook, the 
polynomial dispersion equation (17th degree) has been transferred to MatLab for numerical study. A systematic search 
of the most unstable modes has thus been conducted in terms of a, 75 and 17 s for R — 1/1836 and 1/100. For each 
parameter set (a, 7b, Ob), the program returns the Z z and Z x components of the most unstable mode together with 
the real and imaginary parts of the corresponding complex frequency u> r +18. 

The function 8m = max{(5, (k x , k z ) £ R 2 } is clearly continuous. But the real part u>m corresponding to Sm, together 
with the wave vector Zm defining the mode growing at Sm, may well evolve discontinuously as the most unstable 
mode can perfectly "jump" from one location of the 2D spectrum to another (Bret et al. 2008). For example, the 
Z z component of the fastest growing mode along the beam axis, switches abruptly from Z z = 1 to 1/a when the 
Buneman mode overcomes the Two-Stream one. Indeed, these discontinuities are a way to determine the domains 
where any given mode governs the system: as long as the parameters evolve in such a way that the very same mode 
keeps dominating, Zm and lom are continuous 2 . But when the dominant mode changes, one of these functions will 
suffer some discontinuity. Within the present cold model, the most unstable oblique modes are often "sent" to Z z = 00 
so that an oblique-oblique transition will not necessarily trigger some discontinuity of the perpendicular component. 
But this is neither the case of the parallel component Z z nor of the real part of the frequency lom- Figure 5 has thus 
been elaborated exploiting these properties and is commented in the sequel. Here again, and for better clarity, we 
single out the diluted beam regime a -C 1 from the symmetric case a — 1. The main reason for such a structure is 
that only these two extremes are analytically accessible, whereas the intermediate case definitely requires numerical 
assistance. 

5.1. Diluted beam 



2 The roots of a polynomial are continuous functions of its coefficients. See Uherka & Sergott (1977) for a proof. 
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Fig. 6. — Low a partition of the phase space defined by Eqs. (22,23,24) for 1 <fl B <R~ 1/3 (see Sec. 5.2). 

Gathering the aforementioned data, the frontiers between the domains governed by different modes are sketched 
in terms of (7b, a) on Fig. 5. Although a is limited to the range [0.1,1], the lower part of the graph can be easily- 
deduced from the diluted beam regime expressions reported in Table 1. This process is repeated for various Ob's and 
#= 1/1836,1/100. 

Starting with £Ib =0, the findings of Bret & Dicckmann (2008b) are here confirmed numerically. Filamentation 
dominates for large beam-to-plasma density ratios. The frontier visible in the lower-right corner pertains to the 
Obliquc/Buneman transition. The growth rate of the un-magnetized oblique competing mode in this case 3 reads 
V / 3/2 4 / 3 (a/7b) 1 / 3 (Famberg et al. 1970). The equation of the border is therefore simply, 



a = Rj b . 



(21) 



For a — R and 7$, = 1 , Eqs. (10,11) arc identical so that this frontier extends all the way down to a = R. Bellow this 
line, Buneman modes govern the linear phase. 

From = 1, the evolution is three- fold. To start with, the Two-Stream instability governs an increasing weakly 
relativistic region (lower-left corner) as its magnetized growth rate remains the same, while oblique modes are made 
less unstable. By virtue of the same kind of effects, Buneman modes gain weight on the lower-right corner because 
their growth rate too does not vary with the magnetic field. Finally, Filamentation domain shrinks as the magnetic 
field progressively shuts it down (Cary et al. 1981; Stockem et al. 2007). Regarding the border of the Two-Stream 
region, a look at the Z x component shows it switches from to 00 when crossing the limit. A comparison between the 
growth rate values in this region and the expressions gathered in Table 1, shows that the Two-Stream mode competes 
with the Upper- Hybrid-Like one located at Z z <~ &,B/"fb+ \/l + O s . The Two-Stream/Oblique frontier is thus defined 
for a <C 1 by, 



V3 a 1/3 _ sja 



2 4 / 3 



76 



20, 



<^> a = 



47 b 6 



Regarding the Buneman/Upper-Hybrid-Like frontier, the growth rate 
Buneman mode overcomes it. This transition occurs thus for, 



(22) 

\yfa/£lb simply decreases with a until the 



dI/3 _ 

2 4 / 3 2fl B 



3i? 2/3 o2 



(23) 



Below this intersection, Two-Stream modes are directly competing with the Buneman ones. This later frontier is thus 
defined by 

a = R-yl (24) 

The partition defined by these 3 equations in the diluted regime is thus pictured on Figure 6 for O^ larger than 1 but 
smaller than <~ i? -1 / 3 (see Sec. 5.2 for this limitation). 

For larger values of f2s, the previous trends amplify. The oblique domain (central region) is bounded towards the 
low en's because oblique growth rates are scaled like a 1 / 2 (or a 1 / 3 ), while the Buneman growth rate does not vary with 
a. The evolution is of course faster with R = 1/100 than 1/1836 because the Buneman grows faster in the second 
case. 

From CIb ~ 5.8, the Buneman modes gain more "territory" according to the aforementioned mechanisms. Note- 
worthily, new modes appear to govern some portion of the phase space at rather high density ratio and moderate 
Lorentz factor. For R = 1/1836, the inserted spectrum maps on Fig. 5 show that a small part of the phase space (blue 
one) pertains to the oblique mode described in Section 4.3. This "middle region" is analytically involved to explore and 
such results could only be derived numerically so far. The same mode governs a similar region for R = 1/100, but in 
that case, another kind of oblique mode (green regions) also intervenes. Why can such modes lead the linear phase for 



3 The "Oblique s0 " of Table 1. 
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Fig. 7. — (Color online) Growth rates S (u) p units) of the competing modes in terms of tj,, for a = 1, R = 1/100 and 1/1836 (numerical 
evaluation). The merged Two-Stream/Buneman modes (TSB, in red) varies like 7^ 3 ^ 2 for 75 < R~ 1 / 3 , and like 7J -1 after. Oblique modes 
(in blue) vary like 7 ( J~ 1 ' 15 and Upper-Hybrid-Like modes (UHL, black) like HJ^ 1 . Due to the non-monotonic behavior of Oblique modes for 
moderate 75, the determination of the hierarchy is highly non-trivial within this energy range. 

R = 1/100 and not R = 1/1836? The solution lies in the mode coupling already observed in Sec. 4. Figures 3 (A,C) 
show how the Buneman modes interfere with the rest of the spectrum according to the value of R. For R — 1/100, 
interferences are stronger than for R = 1/1836 because the spectrum is wider. While the beam is diluted (as is the 
case for these plots), coupling effects are limited because the various modes are well separated in the k space. As the 
beam-to-plasma density ratio a increases, the Buneman spectrum progressively merges with the rest of the unstable 
modes. Interferences becomes thus potentially more intense, and all the more that R is large. It turns out that the 
green region in Fig. 5 arises from a mode coupling which is much less excited for R = 1/1836 than for R — 1/100. We 
find here that tuning the electron-to-proton mass ratio in order to speed up the system dynamic may bear qualitative 
consequences by changing the nature of the dominant mode. 

5.2. Diluted beam, ultra-magnetized regime 

The outcome of the ultra-magnetized regime is relevant for Pulsars physics, for example. For diluted beams, an 
increasing magnetic field progressively stabilizes every modes except the Two-Stream, the Buneman and the Oblique 
one (see Table 1). Both Two-Stream and Oblique are scaled like a 1 ^ 3 /^, but the pre-factors favors Two-Stream modes. 
Competition is eventually between Buneman and Two-Stream with a frontier defined by Eq. (24). The magnitude of 
£Ib needed to trigger this picture can be derive by setting a = 1 in Eq. (23) as if the "triple-point" pictured in Fig. 6 
was rejected to a = 1. Such criterion places the ultra-magnetized regime beyond, 

t^-iT 1 / 3 , (25) 

which fits well what is observed on Fig. 5. Let us now detail the symmetric case and detail the mode hierarchy in 
terms of (7^, fis). 

5.3. Symmetric case, a = 1 

When a reaches unity, Table 2 indicates that the competition is between the merged Two-Stream/Buneman, the 
Upper-Hybrid-Like and the Oblique modes. However, a comparison of the growth rates is possible as long as these 
modes are well separated in the k space. When overlapping, the growth rate interfere, and analytical predictions 
become very involved. A closer look shows that the region, 

fl B < V2R lb , (26) 

is thus defined. This limit lies bellow the cancelation threshold for the Filamentation instability, Hb — (3y/2r/b, at 
least up to 7b = 1 /R. Indeed, the numerical exploration of this region shows that Filamentation governs the system 
almost as long as it is not canceled. This allows for the tracing of the Filamentation domain in Figure 8, together with 
the other limits explained in this section. As Filamentation vanishes, it evolves smoothly into the Upper-Hybrid-Like 
modes with ku k Z z / 0. 

When condition (26) is fulfilled, it is possible to directly compare the growth rate of the three competing modes. The 
resulting situation is rendered on Fig. 7 and varies strongly with R. The origin of such i?-sensitivity is the fact that 
only the Two-Stream/Buneman growth rate depends on the mass ratio in the large 75 limit. While its 75 variation is 
monotonous, the variation of the Oblique modes is not, rendering the determination of the hierarchy highly non-trivial 
for moderately relativistic beam. 

Fig. 7 shows that for R = 1/100, Oblique modes are only allowed to overcome the others as long as is not too 
large. Numerically, it is found that for fls > 300, they are overcome by Two-Stream/Buneman regardless of the beam 
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Fig. 8. — Modes hierarchy for a = 1, plotted for R = 1/100 and 1/1836. Filamentation modes (sometimes referred to as "Weibel") govern 
almost as long as their are not canceled. The black star and the square pertain to the parameters chosen for the Intergalactic Streams 
and Relativistic Shocks described in Sees. 6.2 & 6.3 respectively. The marks have been placed on the 7j, axis because the systems they 
represent are un-magnetized. 



represent arc un-magnetized. 

energy. Beyond this threshold, Upper-Hybrid-Like modes govern 11 



n B < 



2 7/6 



lb 



(27) 



The resulting hierarchy map is display on Figure 8 and uncover an intriguing Bubble-like Oblique domain. 

When decreasing R down to its realistic value 1/1836, Fig. 7 shows how the large jb part of the Two- 
Stream/Buneman curve is shifted down by virtue of its R 1 ^ 6 /^ scaling (see Table 2). Even for Clg = 500, the 
Oblique growth rate (plain blue curve) surpasses the Two-Stream/Buneman from Yf, ~ 5 to ~ 200. Whether such 
situation holds for any fls is an open question, but it has been checked that it does at least up to Hb = 10 3 . In order 
to draw the hierarchy map, let us first determined to Lorentz factor for which Two-Stream/Buneman modes overcome 
the Oblique one. The corresponding value of y& fulfills, 



1 1 



: 7i 



1.15 



V3 R 1 ^ 
2 7 /6 



lb 



lb 



0.055 



(28) 



yielding y& ~ 234 for R = 1/1836. The hierarchy map here is the one pictured on Fig. 8 for R = 1/1836. For 
Vis > 2 x 234 1 15 = 1060, the system is governed with y& increasing by Two-Stream/Buneman modes, then Oblique, 
then Two-Stream/Buneman again, then Upper-Hybrid-Like and finally, Filamentation. For smaller CIb's (yet, larger 
than ~ 4), the last Two-Stream/Buneman step is just skipped, as the system goes directly from the Oblique to the 
Upper-Hybrid-Like regime. 

It is therefore uneasy to define here some ultra-magnetized regime. Unlike the diluted beam case where only two 
modes are left to compete for a reasonably high magnetic field defined by Eq. (25), we find here that every modes 
involved are likely to play a role, regardless of the magnetization. For R = 1/100, one could defined a so-called 
ultra-magnetized regime from the top of the Oblique "Bubble" located at Qb ~ 300. For larger magnetization, only 
three modes are left to compete. But such definition is no longer possible for R — 1/1836. 

Finally, let us add a remark regarding the accuracy of the borders calculations. While it is perfectly possible to 
numerically compute hierarchy maps 5 and 8 with higher definition, such progress is not imperative because two 
different systems located on each side of a given border are likely to evolve in a quite similar way during the linear 
phase. Let us assume the border between two mode domains A and B goes through a — oto, in such a way that A 
dominates for a < a®. For a — the growth rates 5 a and 8b are strictly equal. If a = ao ± e, then 5 a — 8 B = T^'- If 
both A and B are excited with the same initial amplitude, the time required for the fastest to overcome the slowest only 
by a factor e is ~ lu^ 1 je' . Depending on the system considered, this time may well exceed the duration of the linear 
phase for s' small enough. It would thus be inappropriate to claim that only A or B are relevant on their respective 
side of the border. Even if the dominant mode evolution can be discontinuous, the system dynamic stemming from 
the growth of the entire spectrum should evolve smoothly when crossing a frontier. One can focus on either A or B 
only far away from the borders so that the precise calculation of its location is not crucial. 

6. APPLICATIONS 

We finally turn to the determination of the modes hierarchy in terms of the parameters a, Yi> and Qb for the three 
settings mentioned earlier. Parameters considered in each case are reported in Table 3 and have been chosen after 
Karlicky (2009), Lazar et al. (2009) and Silva et al. (2003). The respective growth rate maps have been plotted on 
Figures 9(A-C). Although the present theory is limited to the cold case, possible kinetic effects are indicated in each 
case, and even calculated for the Solar Flares environment. 
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TABLE 3 

Parameters used for the Applications. 



Setting 



n b /n. 



Dominant mode 



Solar Flares 1 
Intergalactic 2 
Shocks 3 



1/8 



1.34 

1.005 

10 



0, 1 







Obliquegoi Two-Stream 

Oblique_eo 

Filamcntation 



1 
1 



Trom Karlicky (2009), = 2/3 with R = 1/1836. 
2 From Lazar ct al. (2009), /3 = 1/10 with R = 1/1836. 

3 From Silva et al. (2003). Two colliding pair plasmas without background ions, i.e. R = 0. 



6.1. Solar Flares 



We here perform the calculation for one magnetized and one non-magnetized case. The system location in the 
parameters phase space is designated by the red crosses on Fig. 5, and the growth rate maps are displayed on Figs. 
9(A). For both scenarios, the Buneman modes are clearly visible at Z z <~ 8 = 1/a, but cannot compete with the 
Two-Stream ones at Z z ~ 1 because the beam is not relativistic enough. With such a beam-to-plasma density ratio, 
a much higher relativistic factor is required for the Two-Strcam/Buneman transition. 

For the non-magnetized case (upper Fig. 9A), this weakly diluted beam system is governed by non-magnetized 
"Obliqueso" modes a.t Z z ~ 1. Although not dominant, the Filamentation instability 4 plays an important role as 
its growth rate is quite close to the largest one. The introduction of the magnetic field (lower Fig. 9A) damps 
Filamentation as well as oblique mode, resulting in a Two-Stream driven system. Noteworthily, Figure 5 shows that a 
small variation of the Lorentz factor can trigger a dominant mode transition for Qg = 1. For R = 1/100, the system is 
close to a triple-point where Two-Stream, Buneman and Upper-Hybrid-Likc modes grow the same speed. The system 
behavior in this case is therefore very sensitive to the parameters choice. 

A word of caution is needed here before concluding with respect to the Solar Flares environment. The simulation 
performed in Karlicky (2009) accounts for an electronic temperatures of 21.4 x 10 6 K (21.4 MK). Indeed, analysis of 
soft and hard X-ray flare data indicate that plasma temperatures up to 40 MK are obtained in large flares (Aschwanden 
2002). Furthermore, velocity distributions can be highly anisotropic due to stronger heating along the magnetic field 
lines (Fisk 1976; Miller 1991; Miller et al. 1997). For a mode with wave vector k growing at growth rate Wj, a 
temperature spread Av can be neglected providing IAvIc^ 1 <C k~ x (Famberg et al. 1970). Such condition simply 
ensures that during one e-folding time, particles are traveling almost the same path when compared to the wavelength 
considered. Focussing on the fastest growing Obliqueso mode for the non-magnetized case (upper Fig. 9A), we find, 



Accounting for |Av| = 0.08c (i.e. T = 20 MK), the condition above translates 0.08 <C 0.31. The only change when 
considering the Two-Stream dominated magnetized case (lower Fig. 9A) is that the Lorentz factor on the right-hand 
side has the exponent -1 instead of -1/3. In such case, the condition reads 0.08 <C 0.25. The cold approximation seems 
therefore reasonable here, although condition (29) is only weakly fulfilled. An interesting consequence of Eq. (29) is 
that the condition is not homogenous over the unstable spectrum. A kinetic theory of the full unstable spectrum is 
currently under elaboration, together with a rigorous mode-dependent definition of the kinetic-fluid transition (Bret 
& Gremillet 2009). For a clearer assessment of the accuracy of the cold approximation in the present case, the insert 
in the upper Fig. 9A displays a kinetic calculation of the most relevant portion of the unstable spectrum, assuming 
electronic beam and plasma temperatures of 20 MK. Such calculation has been performed neglecting ion motion, which 
is perfectly valid here. One can check how the portion of the unstable spectrum only slightly departs from its cold 
counterpart, while the largest growth rate switches from 0.3u> p down to 0.22w p . Note that similar calculation for the 
magnetized case is yet to be done. 



The corresponding spectrum is plotted on Fig. 9B. With a small (3 and a beam-to-plasma density ratio of 1, this 
case is also sensitive to the parameters (see the star on Fig. 8) . We find here that non-magnetized "Obliqueso" modes 
govern the system, though Filamentation is far from being shut-down. But such a weakly relativistic symmetric system 
can switch between Oblique or Filamentation regimes trough very small parameters variations. When accounting for 
kinetic effects, Two-Stream/Buneman modes can also compete because they are the least sensitive to temperature. 
Indeed, for R = 0, these three modes grow almost exactly the same way for a beam with ji, = 1.1 and temperature 
100 keV (Bret et al. 2008). Such result was found for a plasma temperature of 5 keV, but similar conclusions can be 
drawn with different plasma thermal spread. The discussion at the end of Section 5 is relevant in this case as this 
system is bordering various frontiers. 

6.3. Relativistic Shocks 

These kind of structures are currently extensively studied by means of PIC simulations (Medvedev & Loeb 1999; 
Silva et al. 2003; Nishikawa et al. 2003; Milosavljevic et al. 2006; Chang et al. 2008; Spitkovsky 2008a,b; Martins et al. 




(29) 



6.2. Intergalactic Streams 



4 Denoted "Weibel" in Karlicky (2009). 
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Fig. 9. — (Color online) 2D unstable spectrum in terms of the reduced wave vector Z for parameters mentioned in Table 3. (A) Solar 
Flares. (B) Intergalactic streams. (C) Relativistic shock. The insert on the upper left figure shows the result of a fully kinetic calculation 
with R = 0, assuming electronic beam and plasma temperatures of 20 MK. 

2009), due to their role in the Fireball model for Gamma Ray Bursts (Piran 2000). The typical scenario arising from 
these studies is two-fold: to start with, two relativistic plasmas shells collide and the resulting instability eventually 
generates a quasi-steady propagating shock. These plasma shells can be pair plasmas (Silva et al. 2003; Chang et al. 
2008), or electron- ion plasmas (Spitkovsky 2008a; Martins et al. 2009). The un-magnetized pair plasma case can be 
analyzed here setting R = £Ib = 0, and the system is found clearly governed by the Filamentation instability (see 
Fig. 9C). This conclusion supports the emphasis put on this instability by the aforementioned authors 5 . Moreover, 
such domination is quite robust: as long as nb/n p = 1, 7b = 10 and Qb — 0, the only alternative to the Filamentation 
instability are the non-magnetized oblique modes. Accounting for kinetic effects, a transition to such a regime demands 
a beam temperature of the order 10 6 keV (Bret et al. 2008). Figure 8 also suggests that some magnetized version of 
the system with Qb > 7&'V / 2i? could trigger a transition to the magnetized oblique regime. But the bigger threat to the 
domination of Filamentation seems to be the beam-to-plasma density ratio. As evidenced by Fig. 5, this instability 
dominates the ultra-relativistic regime only for strictly symmetric systems. At 7b = 10 and 100, Filamentation governs 
only for nb/n p > 0.6 and 0.85 respectively. The symmetric system hypothesis has been so far related to the simulation 
of equivalent density colliding shells, within the Fireball framework. Since a slight change in the density ratio may 
trigger a dominant mode transition, subsequent work should be needed to test the robustness of the shock formation 
in this respect. 

The second stage of the scenario involves particles which are accelerated by the shock (Spitkovsky 2008b), escape 
it, and interact with the upstream medium generating more instabilities. Considering now Q,b = 0, 7b = 100 and 
a = 10~ 2 , Fig. 5 shows that Oblique^n modes dominate for R = 1/1836 while Buneman modes do so for R = 1/100. 
A more detailed evaluation of the spectrum should be done from this point, accounting for the non-thermal energy 
dispersion of the escaping particles, and for the upstream ion temperature. Sorting out this issue may be important 
because an interaction governed by the purely electrostatic Buneman instability is less likely to feed the magnetic 
turbulence needed in the Fireball scenario for synchrotron radiation emission. At any rate, the small beam-to-plasma 
density ratio implied in this second stage should prevent the Filamentation instability from playing any prominent 
role in the linear phase. 

7. DISCUSSION AND CONCLUSION 

An exact cold fluid model for a relativistic beam-plasma system including a guiding magnetic field and ion motion 
has been implemented, allowing for a unified description of every possible instabilities arising within such systems. 
After briefly reminding, or deriving, the key results for each kind of unstable modes, the hierarchy between them 
has been established in terms of the beam-to-plasma density ratio, the beam Lorentz factor and the magnetic field 
strength, considering electron-to-proton mass ratios of 1/1836 and 1/100. 

In the diluted beam regime, the competition is mainly between the Two-Stream, the Buneman and the Upper- 
Hybrid-Like modes, as the hierarchy diagram in terms or (7^, a) adopts a typical "V" shape for a given Qb (though 
not too high) . Two-Stream and Buneman take advantage of the magnetic field which leave both of them unaffected 
while stabilizing the rest of the spectrum. Additionally, the Buneman modes, arising from the interaction of the 
electronic return-current with the background protons, govern the highly relativistic regime where the Two-Stream 
modes are strongly reduced. The ultra-magnetized regime with > i? -1 / 3 is eventually governed by either the 

5 Here again, "Filamentation" is generally labeled "Weibel" in the literature. 
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Two-Stream or the Buneman instabilities. 

When leaving the diluted beam region, the shape of the hierarchy diagram is far from being trivial, and two different 
kinds of Z z ^ modes may dominate. The first ones are the ones already governing a part of the parameters space 
from VLb = 0. In the cold fluid un- magnetized model, they are found at Z z ~ 1 and Z x — oo, whereas temperature 
effects would give them a finite perpendicular component. Under the action of the magnetic field, they continuously 
evolve to the Upper- Hybrid-Like modes refereed to in Table 1, with growth rate oc y/a/fls- But from fls ~ 4.3 and 
3.6 for R = 1/1836 and 1/100 respectively (numerical evaluations), another kind of oblique mode is likely to govern 
the spectrum. Unlike the previous ones, these modes reach their maximum growth rate for one single wave vector and 
have the growth rate scaled like a 1 / 3 With such a 7b scaling, their domain shrinks with increasing beam energy, 
as observed on Figs. 5. 

The transition is thus made with the symmetric regime a = 1. Here, Two-Stream and Buneman modes merge 
and dominate at low Lorentz factor. According to Fig. 8, there is no such thing as a simple description of the 
ultra-magnetized regime. For a large mass ratio R = 1/100 and > 300, Two-Stream, Upper-Hybrid-Like and 
Filamentation modes are likely to dominate, depending of the Lorentz factor. At lower magnetization, Oblique modes 
are likely to play a role for moderate beam energy. At lower mass ratio with R = 1/1836, the hierarchy is even more 
complex as the Oblique "Bubble" found at higher R expands towards the large Qb regime. 

While the ultra-magnetized regime is governed by electrostatic modes with flow aligned wave- vector, this is not true 
for density ratios close to one. Indeed, the later can be determined by virtually any kind of modes. In this respect, 
diluted systems are simpler to analyze, or simulate, than symmetric ones. 

Some modes, like the Buneman's, strongly depends on the electron-to-proton mass ratio R while others don't. As a 
result, the hierarchy is i?-dependent, and lowering the ratio for computational purposes can bring some qualitative, not 
just quantitative, changes. Figure 5 shows for example that a system with a = 0.1, 7b = 20 and il^ = 1 is governed 
by Oblique modes for R = 1/1836 and Buneman ones for R = 1/100. Equation (24) shows that the nature of the 
diluted beam ultra-magnetized regime strongly depends on R. 

By virtue of its possible discontinuous nature, as explained in Sec. 5, a change of the dominant mode can have 
dramatic consequences. One of them can be illustrated through the Two-Stream/Buneman transition defined by Eq. 
(24) for diluted beams. Because Two-Stream modes are found near Z z = 1 and Buneman ones near Z z = l/a, the 
typical size of the structures generated is co p /vb on the left hand side of the border, and (n p /nb)u> p /vb on its right side. 
As long as the beam is not too diluted, these two quantities remain similar. But for a 10 3 times diluted beam, both 
kind of structures differ by 3 orders of magnitude. For a relativistic symmetric system with a = 1 and 7b > ^i? -1 / 3 , 
a transition from the Two-Stream/Buneman regime to the Upper-Hybrid-Like one implies a switch from a Z z = \/2R 
dominant mode to a Z z = Qb/% one which again can span several orders of magnitude. 

These transitions have to do with the parallel wave-vector component rather than with the perpendicular one. 
Similar transitions have been recorded in this paper for this later component, but unstable modes located at Z x = oo 
acquire a finite normal component as soon as temperature are accounted for. While kinetic effects may temper the Z x 
transitions, such should not be the case for the magnitude of the Z z transitions, because parallel components are less 
sensitive to temperature that normal ones. 

Turning now to the astrophysical settings considered here, it is worth stressing that the value of the present theory 
consists more in pointing what the main instabilities could be, rather that confirming what they actually are (at least 
within the current models). Considering Solar Flares of Relativistic Shocks for example, it has already been checked 
that the dominant instabilities are the ones stressed in Sec. 6. Even in Filamentation, or Weibel, modes are sometimes 
designated instead of Oblique ones, these latters are increasingly specifically discussed in the literature, and identified 
as such (Jaroschek et al. 2005; Niemiec et al. 2008; Ohira & Takahara 2008; Frederiksen & Dieckmann 2008; Kong 
et al. 2009). But what PIC simulations can hardly do is indicating extensively which kind of modes could take the 
lead when changing the parameters, including the electron-to-proton mass ratio. 

For Solar Flares, the domination of Obliqueeo modes is quite robust in the non-magnetized case. Although the 
cold approximation should correctly describe the most unstable mode, given the temperatures involved, a fully kinetic 
calculation of the spectrum has been done to confirm this point (see upper Fig. 9A). The situation is quite different 
when magnetizing the system. For the magnetic parameter considered here, the lower Fig. 9A clearly show how 
some Uppcr-Hybrid-Like modes at Z z ~ 2 can also play a role. Indeed, Fig. 5 shows that for a realistic value a the 
electron-to-proton mass ratio R — 1/1836, the transition only needs a slightly higher relativistic factor. In such a 
case, the instability will heat the system both in the parallel and perpendicular direction while the Two-Stream driven 
counterpart preferentially heats it along the beam direction (Karlicky 2009). Furthermore, tuning the mass ratio to 
speed up the dynamic can here bring qualitative consequences, as evidenced in Fig. 5. It is thus found that the 
evolution of a magnetized electron beam in this setting, together with the kind of heating provided, is quite sensitive 
to the parameters. 

Regarding Intergalactic Streams, the point made in Lazar et al. (2009) is that given the parameters involved, 
Filamentation instability grows faster than the rest of the spectrum. The expected saturation level of the magnetic 
field thus produced is then found consistent with the measurements. In this respect, Fig. 8 shows that the cold version 
of the system is bordering the Two-Stream/Buneman-Filamentation frontier, so that Filamentation "leadership" is 
eventually not so clear. Admittedly, kinetic effects are important as Filamentation is found to dominate through its 
interaction with some temperature anisotropy, but such growth rate enhancement has been found to operate in other 
parts of the spectrum as well (Bret et al. 2005). The question remains open as to know how strong Filamentation 



Weibel, Two-Stream, Filamentation, Oblique, Bell, Buneman. 
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leadership needs to be for the linear phase to saturate accordingly. Such criterion would likely define some parameters 
window in which Weibel- like instabilities could be responsible for intergalactic magnetic fields. 

Finally, the instability-based scenario for shocks formation within the Gamma Ray Bursts framework has been 
examined (Fireball model). Beam-plasma instabilities have been found so far to play a key role at two levels. First, 
they seem to prompt the shock formation itself, as the collisionless encounter of two plasma shells is unstable. Second, 
once the shock has been formed, it accelerates particles through diffusive Fermi acceleration. Particles escaping the 
shock upstream interact with the medium, and generate instabilities responsible for the magnetic turbulence needed 
to trigger synchronic radiation emissions. The first unstable system has a beam-to-plasma density ratio close to unity. 
Indeed, it is necessarily unity in the recent PIC simulations where a beam is reflected against a wall, and eventually 
interact with itself. For such settings, Filamentation instability has been found leading the linear phase, and the 
shock formation has been so far analyze accordingly (Medvedev & Loeb 1999; Spitkovsky 2008a, b). But the present 
calculations suggest that a density ratio slightly smaller than 1 may result in a quite different dominant instability. It 
would be very interesting to test the robustness of the shock formation scenario in this respect. Once the shock has 
been formed, instabilities upstream are thus expected to generate magnetic turbulence. We now deal with a highly 
relativistic diluted beam/plasma system, which should definitely not be governed by the Filamentation instability (see 
Fig. 5). Kinetic effects are unlikely to modify this picture, as Filamentation modes are usually their first "victims" 
(Bret & Grcmillct 2009). Here again, a parameter window is defined allowing for the development of the best magnetic 
turbulence generators: Figure 5 clearly shows how Obliquesoj Upper-Hybrid-Like, Two-Stream or Buneman modes 
arc likely to shape the linear phase on this second phase. Among these four candidates, the last two are reputed purely 
electrostatic instabilities, which should thus be avoided. An accurate characterization of the beam escaping the shock 
upstream will thus be needed to assess the validity, or set limits, to this scenario. 
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APPENDIX 
TENSOR ELEMENTS 
The tensor mentioned in Eq. (6) can be cast under the from, 
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